Specific Heat of the Dilute Ising Magnet LiHo^Yi-^^ 
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We present specific heat data on three samples of the dilute Ising magnet LiHo a: Yi_ :r F4 with 
x — 0.018, 0.045 and 0.080. Previous measurements of the ac susceptibility of an x = 0.045 sample 
showed the Ho 3+ moments to remain dynamic down to very low temperatures 0, 0] and the specific 
heat was found to have unusually sharp features P, 0|. In contrast, our measurements do not 
exhibit these sharp features in the specific heat and instead show a broad feature, for all three 
samples studied, which is qualitatively consistent with a spin glass state. Integrating C/T, however, 
reveals an increase in residual entropy with lower Ho concentration, consistent with recent Monte 
Carlo simulations showing a lack of spin glass transition for low x Q. 
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Extensive work has previously been done to under- 
stand the spin glass transition found in disordered mag- 
netic systems [5J and changes in behavior as the concen- 
tration of magnetic moments is reduced jf|. The mate- 
rial LiHoj;Yi_ x F4 is a nearly perfect example of a dilute, 
dipolar-coupled Ising magnet and is therefore an ideal 
system for experimentally testing theories of simple, in- 
teracting spin models. Despite the apparent simplicity of 
this system's underlying model, however, a series of sur- 
prising results and fascinating effects has emerged from 
the material's rich phase diagram, especially at low con- 
centrations of magnetic Ho 3+ ions 1 

At x = 1 the system has been found to order ferro- 
magnetically with T c — 1.53 K [8|, but below a certain 
amount of dilution (x ~ 0.25), there is enough random- 
ness and frustration (due to the angle-dependent dipolar 
interaction) that the system becomes a spin glass jl|, Q . 
At x = 0.045, however, ac susceptibility experiments 
have shown the material not to freeze down to very low 
temperatures [H- The absorption spectrum x"( w ) 1S ob- 
served to narrow with lower temperature where typically, 
in a spin glass, the absorption spectrum becomes wider 
as freezing of the moments leads to longer relaxation 
times 10]. Furthermore, at temperatures below 100 mK, 
there appears to be a gap in the absorption spectrum and 
coherent, low frequency oscillations with lifetimes of up 
to 10 s are observed. These effects have been attributed 
to clusters of roughly 260 Ho ions acting as largely inde- 
pendent oscillators [2j. 

Since the dipolar coupling between the Ho moments is 
a long-range interaction, it has long been theoretically ex- 
pected that there should be no finite concentration of mo- 
ments (x) at which the ordering (or freezing) temperature 
of the system drops to zero [ll| . The ac susceptibility ex- 
periments performed on the x = 0.045 sample, however, 
seem to contradict this theory as there is no sign of freez- 
ing even down to 50 mK. Recent Monte Carlo simulations 
of dipolar-coupled Ising moments randomly placed on a 
cubic lattice do suggest that there is no spin glass transi- 
tion for x < xc — 0.20 which could explain the unusual 
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FIG. 1: Total measured specific heat for x = 0.018, x = 0.045 
and x = 0.080. The solid line is the total non-interacting 
specific heat calculated by diagonalizing the crystal field and 
nuclear hyperfine Hamiltonians. 



dynamics seen at x = 0.045 |4|. There is also evidence 
in recent /xSR experiments on LiHoo.o4sYo.955F4 of spin 
dynamics persisting down to low temperatures [l2l ]. 

This unique spin liquid or "anti-glass" state has also 
exhibited unusually sharp features in its specific heat at 
around 110 mK and 300 mK. These features were qualita- 
tively reproduced in numerical simulations using a model 
based on quantum entanglement of pairs of moments and 
a pair-wise 'decimation' procedure in which the sharp 
heat capacity features correspond to maxima in the dis- 
tribution of dipolar couplings in the system This 
simulation was also able to reproduce a T~ a75 behavior 
of the dc susceptibility which was observed experimen- 
tally. It is not known whether there is a relation between 
these heat capacity signatures and the anomalous dy- 
namics observed by ac susceptibility. 

In this letter we present specific heat data taken on 
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three stoichiometrics in this series: x = 0.018, 0.045 and 
0.080. The samples studied are high quality single crys- 
tals grown with the Bridgman technique [29| • Crystalline 
quality was verified by high resolution diffraction on a 
fine focus Cu rotating anode generator equipped with 
a high resolution Ge (220) four-crystal monochromator 
and a Huber 4-circle diffractometer. The measurements 
revealed extremely sharp Bragg peaks (6*fwhm < 0.015°) 
for all reflections, indicating high crystalline perfection. 
No twinning was observed. Extensive diffuse scatter- 
ing measurements revealed no diffuse scattering near or 
away from the Bragg peaks, or satellite peaks that could 
be associated with any disorder or short range ordering. 
Small ~ 100 /im fragments were taken from each sam- 
ple and crystallography data sets were measured using 
a molybdenum rotating anode, kappa diffractometer and 
CCD area detector. All three data sets refined well with 
Ho substituting for Y in the expected tetragonal (74i / a) 
structure 

Measurements were performed using the quasi- 
adiabatic method with a long time-constant r of relax- 
ation. No substrate was used in these experiments and 
the heater, thermometer and weak link were glued di- 
rectly onto the sample which was suspended from very 
fine nylon threads. Using a substrate can lead to an un- 
derestimate of the heat capacity due to large thermal 
resistances between the sample and substrate. The ad- 
dendum was determined to be less than 0.1% of the sam- 
ple's heat capacity. Samples were typically discs ~ 8 mm 
in diameter and ~ 1 mm thick. 

A Ru02 resistor (1 kf2 at 300 K) was used as a ther- 
mometer and a 10 kf2 metal-film resistor was used as a 
heater, both with thinned alumina substrates. Leads to 
the thermometer and heater were 6 /im diameter, ~ 5 
mm, NbTi, superconducting wires with a thermal con- 
ductance of i^NbTi — 8 x 10 -11 W/K at 1 K and at 
least a factor of 10 smaller at 100 mK. The thermal con- 
ductance from the thermometer and heater to the sample 
(Kts and Kyis respectively) were measured to be greater 
than 10~ 8 W/K at very low T (< 50 mK). The weak link 
connecting the sample to the dilution refrigerator mixing 
chamber was made from manganin wire and had a ther- 
mal conductance -Kwl — 1 x 10~ 7 W/K at 100 mK. Cal- 
culations show that the temperature of the thermometer 
differs from that of the sample by less than 0.1%. 

Thermometer resistance measurements were made 
with a LR-700 AC Resistance Bridge. The cell was con- 
tained in a copper radiation shield and the cryostat was 
surrounded by a lead shield and two /i-metal shields to 
attenuate the external magnetic field. The thermometer 
resistance was consistent with a standard RuC>2 temper- 
ature dependence 14| with no indication of self-heating 
in the range of our data. The RuC>2 thermometer was 
calibrated to a calibrated LakeShore Ge resistance ther- 
mometer and a CMN susceptibility thermometer. 

Time constants on the order of several hours (at the 
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FIG. 2: Diagram of relevant thermal links in experimental 
apparatus (a) and an example heat pulse showing linear fits 
and extrapolation to the mid-point of the pulse (b). 



lower temperatures) ensured that the sample was cooled 
very slowly and was therefore able to reach equilibrium. 
Cooling the sample even more slowly did not have a no- 
ticeable effect on the measured heat capacity. Temper- 
ature data was collected for up to 30 minutes on cither 
side of the heat pulse and the heat capacity is given by 
C = Q/AT where AT is obtained through extrapolation 
to the midpoint of the pulse as is shown in Figure HJb). 

For temperatures below 1 K, the specific heat of 
LiHo 3; Yi_ 2 ;F4 is dominated by a broad feature which 
arises from the I = 7/2 nuclear-spin degrees of free- 
dom. The single-ion Hamiltonian (neglecting the dipo- 
lar and nearest-neighbor exchange interactions) is given 
by H = Hcf + Hhf + Hq. The 4/-electrons in Ho 3+ 
are tightly bound, resulting in a significant nuclear- 
hyperfine interaction: 7Yhf = AI ■ J. 7Yq is the nuclear 
quadrupole interaction and Hcf — J2i m ^T^T * s the 
crystal field potential caused by surrounding ions (the 
0™'s are Steven's operator equivalents). 

If the crystal field Hamiltonian is diagonalizcd by itself, 
one obtains a ground-state, Ising doublet with an effec- 
tive (/-factor of 13.8 and a next excited state at around 
11 K. It is then, in some cases, safe to assume that the 
electronic moments are perfect Ising spins and the spe- 
cific heat can be expressed as the sum of an electronic 
contribution AC and a nuclear contribution 
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where x m = —Ajg gm/2gjT + Pm 2 /T and m = 
—7/2,. ..,7/2. A fit of this form was successfully ap- 
plied by Mennenga et al. to the specific heat of LIH0F4 
below the transition temperature [HI]. We have made 
corrections to this form by diagonalizing the entire non- 
interacting Hamiltonian (a 136 x 136 matrix). We have 
used Aj/ks — 40.21 mK, determined by EPR exper- 
iments on LiHoo.02Yo.98F4 16]. This is similar to the 
value Aj/ks — 39.8 mK found for the pure mate- 
rial [13] • We have assumed an axially symmetric nuclear 
quadrupole interaction of strength P = 1.7 mK which 
was determined with EPR on free Ho 3+ ions [18| . The 
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FIG. 3: Electronic moments' contribution to the specific heat 
(solid line from Figure Q] subtracted) for x = 0.018, x = 0.045 
and x = 0.08 from this work (filled symbols). Also x = 0.045 
from Ghosh et al. 3] and x = 0.167 from Reich et al. [l[ (open 
symbols) . The solid lines are fits of the form of Equation 



crystal field parameters Bf 1 were taken from [19j. The 
resulting single-ion specific heat is shown as the solid 
line in Fig. Q] and has been subtracted from the data to 
give AC in Fig. [3] This more detailed calculation of the 
non-interacting specific heat is lower than Equation [1] by 
~ 4% near the highest point of the curve. 

As mentioned earlier, the validity of this subtraction 
depends on the assumption that the Ho 3+ ions are per- 
fect Ising moments. This is not entirely the case, as the 
nuclear hyperfine interaction introduces mixing with the 
next excited states. Indeed, the nuclear hyperfine inter- 
action has been shown to strongly effect the magnetic 
ordering of the pure material in transverse field [19( and 
the diluted material in the spin glass regime [2(| . Never- 
theless, the total specific heat should approach this non- 
interacting specific heat at higher temperatures (close to 
1 K). A small phonon contribution to the specific heat 
(oc T 3 ) is also subtracted, estimated from the heat ca- 
pacity of the pure material above 1 K [lj|. 

The crystal supplier provided us with nominally 2%, 
4.5% and 8% concentrations of holmium. Assuming the 
correct nuclear hyperfine component, however, the 2% 
sample appears to be closer to 1.8% as the remaining 
term AC should behave as T~ 2 at higher temperatures 
(close to 1 K). In this way, the 4.5% and 8% samples were 
confirmed to be the correct stoichiometry. 

In all three samples, the remaining specific heat con- 
tribution AC is qualitatively similar and found to be a 
broad feature which is somewhat consistent with the heat 



TABLE I: Fitting parameters for AC for x = 0.018, x = 
0.045 and x = 0.08 from this work and data taken from Reich 
et al. (x — 0.167) The peak temperatures T pca k, the 
relative width of the specific heat curve FWHM/T max and the 
measured residual entropy So (assuming a linear temperature 
dependence at low T) are also given. 



Parameter 


1.8% 


4.5% 


8.0% 


16.7% 


Ex/ha (K) 


0.26 


0.32 


0.29 


0.46 


n 


0.85 


1.43 


0.86 


1.89 


Tpeak (K) 


0.11 


0.13 


0.12 


0.17 


FWHM/T pcak 


1.7 


1.6 


1.7 


1.5 


S /R 


0.31 


0.21 


0.00 


0.18 



capacity of a spin glass [2JJ, |22| • The specific heat of a 
spin glass is not expected to show a remarkable feature at 
the spin glass freezing transition as the critical exponent 
a is often negative, in the range -2 to -4 [23|]. Instead 
of probing the actual freezing transition, the spin glass 
heat capacity is more indicative of excitations above the 
transition. The simplest situation is one excited energy 
state Ei above the ground state having a degeneracy n 
with respect to the ground state degeneracy. We can 
then apply fits of the form 



AC = C 



n(Exlk B T) 2 e- E ^ k * T 



(2) 



to the data. The resulting fitting parameters for this data 
and an x = 0.167 sample measured by Reich et al. [l[ are 
shown in Table |U Clearly the size of the specific heat 
features decreases with decreasing concentration x. The 
peak temperature of the curve, however, is very close in 
all three samples and does not appear to scale with the 
Curie temperature (xTq). 

Numerically integrating AC /T with respect to T gives 
the total amount of entropy released over the temper- 
ature range of these measurements. The total high- 
temperature entropy of an Ising magnet is i?ln2 but 
there may be a residual ground state entropy seen in 
doing this integral. Lower temperature data is required 
in order to confidently observe all the release of entropy 
in the system. We have extrapolated the data to K, as- 
suming a linear temperature dependence, before integrat- 
ing AC/T. Measurements at lower temperatures must 
be made in order to determine the temperature behav- 
ior of the specific heat below the maximum, but many 
past measurements have observed a linear temperature 
dependence in spin glasses 2l|, |22| as described by the 



two-level system (TLS) argument [24 

For the 8% sample, this integral reveals approximately 
all of the total expected entropy. In the 1.8% and 4.5% 
samples, however, our measurement observes a smaller 
percentage of i?ln2: 56% and 70% respectively, leaving 
a significant residual entropy So . The residual entropy for 
each sample is also shown in Table HI These values may 
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be compared to previous measurements for x = 0.167 
where 75% of the entropy was measured (also assuming 
a linear temperature dependence below the peak) and 
x = 0.045 where only 15% of i?ln2 was observed over 
the range of the measurement P, Q • In the case of the 
4.5% and 16.7% samples, So is quite close to the value of 
0.199-R predicted by the SK model of a spin glass (25| . 

The Monte Carlo simulations of Snider et al. on a 
dilute dipolar-coupled Ising system predict residual en- 
tropy at x = 0.20 [H . For x < 0.20 they see no spin 
glass ordering and an increasing So with decreasing x as 
a larger number of degenerate grounds states are avail- 
able. This is the same trend observed in our data, though 
experimentally, the magnitude of the residual entropy is 
much larger. In the real system, the point at which spin 
glass ordering ceases must be lower than x = 0.167 as 
this has been observed to be a spin glass [l[ and may be 
closer to x = 0.080 at which point we have observed the 
release of nearly all the expected entropy. 

The relative broadness of the observed features may 
be parametrized by the full width at half maximum 
(FWHM) divided by the peak temperature T max . The 
three samples studied here give values around 1.7 (see Ta- 
ble H}. This parameter is approximately 1.2 in AuFe [26| 
and 1.5 in Eu x Sri_ x S for example [22j. Typically, the 
maximum in the specific heat of spin glasses is found to 
be approximately 20% higher than the spin glass transi- 
tion temperature which is determined by ac susceptibility 
experiments [B| . If this rule of thumb were to apply here, 
it would give spin glass transition temperatures of 90 to 
100 mK for these three samples. 

We have measured the specific heat of three samples 
at and around a concentration of 4.5% holmium, and 
our measurements do not exhibit the sharp features that 
were seen previously 0, The data sets agree well 
until ~ 300 mK which indicates that there is no error 
in stoichiometry. Below this point, however, there is a 
significant discrepancy. Our data, therefore, also does 
not support the theoretical model presented by Ghosh et 
al. 0] • We have taken great care to rule out any experi- 
mental errors such as decoupling of the thermometer from 
the sample. Recent thermal conductivity measurements 
of LiHoo.04Yo.96F4 also do not show any remarkable fea- 
tures at 110 mK and 300 mK 27]. In some systems which 
have x ~ T~ a the specific heat also shows a simple power 
law behavior with a related exponent [281 ] . It would be in- 
teresting to measure C (T) to lower temperatures to look 
for such an effect. 

This specific heat data is consistent with a spin glass 
in that the observed feature is a broad maximum with no 
pronounced anomalies. However, it is not clear how the 
unusual spin liquid or "anti-glass" state observed with 
ac susceptibility [l|, y] should manifest itself in specific 
heat measurements. Based on the numerical simulations 
of Snider et al. 0] , the measured increase in entropy may 



indicate that the system is no longer a spin glass below 
x ~ 0.08 and is instead a spin liquid with many accessible 
nearly-degenerate ground states. 
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